Solvent Mediated Assembly of Nanoparticles Confined in Mesoporous Alumina 
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The controlled self-assembly of thiol stabilized gold nanocrystals in a mediating solvent and confined 
within mesoporous alumina was probed m-situ with small angle x-ray scattering. The evolution of 
the self-assembly process was controlled reversibly via regulated changes in the amount of solvent 
condensed from an under-saturated vapor. Analysis indicated that the nanoparticles self-assembled 
into cylindrical monolayers within the porous template. Nanoparticle nearest- neighbor separation 
within the monolayer increased and the ordering decreased with the controlled addition of solvent. 
The process was reversible with the removal of solvent. Isotropic clusters of nanoparticles were 
also observed to form temporarily during desorption of the liquid solvent and disappeared upon 
complete removal of liquid. Measurements of the absorption and desorption of the solvent showed 
strong hysteresis upon thermal cycling. In addition, the capillary filling transition for the solvent in 
the nanoparticle-doped pores was shifted to larger chemical potential, relative to the liquid/vapor 
coexistence, by a factor of four as compared to the expected value for the same system without 
nanoparticles. 

PACS numbers: 61.46.Df, 68.08. Be, 61.10.Eq 
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I. INTRODUCTION 

Structures of self-assembled nanoparticles have been 
studied extensively in recent years for their unique cat- 
alytic, electronic, and optical properties. The geometry 
of self-assembled structures can vary from 2D sheets of 
nanoparticles, to spherical monolayers or ID nanowire 
arrangementsi2i^i^ Reduced dimensionality systems, 
in particular nanowires, have been of recent interest 
for their potential electronic and optical propertiesA^ 
There has been additional interest in quasi-one dimen- 
sional structures, such as the arrangement of nanopar- 
ticles into cylindrical monolayers, for their similarity 
with biological systems such as some virus protein coat- 
ings and cell microtubules.^'^ Understanding the pro- 
cess of nanoparticle self-assembly in the different geome- 
tries can be vital to production of functional nanoscale 
structures. Recently there have been several experimen- 
tal and theoretical studies that examine the nanopar- 
ticle self-assembly on flat surfaces via the evaporation 
of macroscopic dropletsiSiSiiSiiiiiSii^ One dimensional 
systems typically require templates, thus bulk methods, 
such as droplet evaporation, are not possible. Recent 
studies of quasi- ID nanoparticle structures within cylin- 
drical templates^i^4i^4i^4iSiiLi& primarily probed the end 
product and did not explicitly measure the evolution of 
the self-assembly process. In order to investigate the self 
assembly process itself within confined geometries, it is 



necessary to perform in-situ experiments with precise 
control of the solvent amount, analogous to controlled 
evaporation of macroscopic droplets of nanoparticles on 
flat surfaces. 

Here we describe in-situ small angle x-ray scatter- 
ing (SAXS) experiments that probed the solvent me- 
diated self-assembly of Au-core, colloidal nanoparticles 
within nanoporous alumina. Experiments were carried 
out within an environmental chamber which allowed pre- 
cise control over the amount of solvent condensed from 
vapor into the porous system. Analysis determined that 
the nanoparticles self-assembled into a cylindrical mono- 
layer that evolved with the addition and removal of liq- 
uid. The evolution of this cylindrical monolayer struc- 
ture was completely reversible upon removal of liquid, 
albeit with strong hysteresis typical in capillary systems. 
In addition to the cylindrical monolayer, isotropic clus- 
ters temporarily formed during desorption of the liquid. 
These clusters disappeared upon complete removal of the 
liquid from the pores. 



II. EXPERIMENTAL 

Anodized alumina membranes were electrochemically 
prepared using a two-step anodization technique de- 
scribed elsewhere liSiSSiSi The aluminum backing layer 
was then dissolved in HgCl2. Pores were opened by 
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FIG. 1: (LEFT) SEM image shows well ordered hexagonal 
packing of alumina nanopores prepared under similar condi- 
tions to samples used for x-ray experiments. (RIGHT) Bright- 
field TEM image of the nanoparticles (arrow) on the wall 
(faint diagonal lines) of alumina pores. The nanoparticles are 
Au-core with octane-thiol coating. Cross-sections were pre- 
pared by ultra micro-tome method from the sample used for 
x-ray experiments. 



floating on phosphoric acid followed by thorough rinsing 
in de-ionized water. The resultant alumina membrane 
consisted of a dense matrix of pores running completely 
through the membrane and perpendicular to its surface. 
The long axes of the pores were parallel and arranged 
in a 2D powder-like arrangement with local hexagonal 
packing. Nearest neighbor distances (center to center) 
of 63±2 nm were determined via x-ray scattering. This 
is consistent with electron microscopy of similarly pre- 
pared samples (see Fig. Alumina nanopores were 
29±4 nm in diameter (from TEM). The macroscopic 
dimensions of the nanoporous membrane were about 
1 cm X 1 cm X 90 microns. 

Nanoporous alumina membranes were further cleaned 
in solvents to remove organic impurities, with 15 minute 
ultra-sonic baths in each of the following (in order): 
HPLC grade chloroform, 99.7% pure acetone, HPLC 
grade methanol, and HPLC grade toluene. The mem- 
brane was then allowed to soak in HPLC grade toluene 
for 24 hours to remove any further impurities. After 
soaking, the membrane was transferred without drying 
to 2 ml of a 2.26 mg/ml solution of octane-thiol (OT) sta- 
bilized Au-core nanocrystal g^^'^"^ in HPLC grade toluene 
for two weeks at 25°C. Fits of SAXS data from bulk 
scattering of a dilute suspension of the same solution 
of nanoparticles used in doping the nanopores yielded 
Rs — 1.2 nm and A/? = 0.16 nm giving a polydisper- 
sity of A]i/Rs — 13%. Since the electron density of the 
Au core is much higher than the organic ligands, Rs was 
primarily a measure of the core radius. 

The nanoparticle solution was drawn into the alumina 
pores via capillary forces^. The initial concentration 
(and volume) of the nanoparticle solution was chosen 
such that the total number of particles in the initial so- 
lution was about 2-3 times the number of nanoparticles 
needed to form a complete monolayer of nanoparticles 
along the walls of the pores (assuming hexagonal close 




FIG. 2: Schematic of the environmental chamber used for 
in-situ x-ray experiments. The outer (inner) cell was kept at 
28° C (32° C) during the experiment. The sample temperature 
Ta was controlled independently from the inner cell and liquid 
solvent reservoir temperature, Tr- The sample was mounted 
in a transmission geometry where the long axis of the pores 
can be rotated (Os) in the horizontal plane. The AT = Ts~Tr 
was controlled to vary liquid condensation on the sample via 
the vapor. 



packing) , but less than that required for complete volume 
filling of the pores. In addition to soaking, the solution 
with the immersed membrane was placed in an ultra- 
sonic bath for 15 minutes per day to facilitate movement 
of the particles into the membrane. The long time pe- 
riod and frequent ultra-sonic baths were imposed to en- 
sure maximal integration of nanoparticles into the porous 
membranes. 

After 24 hours in solution the alumina membrane color 
changed from clear to dark red color, similar to the so- 
lution color. This was an indication that nanoparticles 
had been absorbed into the membrane. The presence 
of nanoparticles in the pores was confirmed with cross- 
sectional TEM of the sample (see Fig. No further 
color changes of the membrane were observed. Addi- 
tionally, no significant color change of the solution was 
observed, indicating that some particles were left in the 
solution. The sample was then removed from the solution 
and allowed to dry in air for five minutes before being 
loaded into an environmental chamberiSiS^ (see Fig. ^ 
for the SAXS experiment. 

Prior to sample mounting, all components of the in- 
ner cell of the environmental chamber were cleaned with 
solvents: 15 minute ultra-sonic baths of each of HPLC 
grade chloroform, 99.7% pure acetone, and HPLC grade 
methanol. During the x-ray experiment the inner cell 
was hermetically sealed and the temperatures of the in- 
ner and outer cells were kept at Tin = 32.0 ± 0.5° C 
and Tout = 28.0 ± 0.5° C, respectively. The actual 
temperature stability during each measurement was bet- 
ter than ±50 mK, though the temperature of the cells 
did vary ±0.5 K between different measurements due to 
the differing heat loads on the sample as described be- 
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low (see the description of AT). The sample mounting 
block was thermally isolated to a high degree from the 
walls of the inner cell with independent thermal control 
via a heater on the mounting block. All temperatures 
were continuously monitored with YSI #45008 thermis- 
tors. A liquid solvent reservoir was created by inject- 
ing 10 ml of HPLC grade liquid toluene into the inner 
cell. Solvent condensation in the sample was precisely 
controlled via the (positive) offset between the sample 
temperature, T^, and the liquid solvent reservoir temper- 
ature, Tr-. AT = T,- T,M^. The chemical potential 
offset (from the liquid/vapor coexistence) is given by: 
A/i « HyapAT/Tr. Here, Hyap = 38.06 kJ/mol is the 
heat of vaporization of toluene^^ and T,- was kept con- 
stant throughout the experiment (as described above) via 
temperature control of the inner cell. Liquid was injected 
into the cell at a large AT w 30 K to avoid rapid con- 
densation in the pores. Since AT could be varied from 
about 50 mK up to 30 K, the corresponding A/i could be 
varied about four orders of magnitude. As AT and thus 
A/i decreased the amount of condensed solvent increased. 
Likewise, solvent was removed by increasing AT. In this 
experiment, the pores became saturated with liquid for 
all AT < 4 K, thus the experimental range of AT probed 
was only 0.5 K< AT < 30 K. X-rays were allowed to pass 
through the chamber via 0.02 mm Kapton© (Dupont) 
windows on the outer chamber and 0.005 inch thick beryl- 
lium windows on the inner chamber. 



III. SMALL ANGLE X-RAY SCATTERING 

In-situ x-ray measurements were carried out at the 
SAXS facility at ChemMatCARS beamline at the Ad- 
vanced Photon Source (Argonne National Lab, Argonne, 
IL). The incident x-ray energy of 11.55 keV was well be- 
low the L3 absorption edge of Au (11.92 keV) to avoid 
fluorescence. The sample chamber was mounted on a go- 
niometer in a transmission geometry with two degrees of 
rotational freedom - sample rotation and tilt - plus the 
standard three translations. The sample rotation was 
done internally and allowed rotation angles, 9^, of the 
sample with respect to the incident beam in excess of 
±90° (see Fig. 1^. A fixed position geometry CCD de- 
tector measured the SAXS intensity at a camera length, 
L, of 1880 mm down-stream from the sample. Powder 
diffraction rings were seen at small angles due to the 2D 
hexagonal packing of the nanopores. Sample alignment 
was performed with the membrane nearly perpendicular 
to the incident beam and thus the pores were parallel to 
the incident beam. The sample rotation and tilt were 
adjusted to maximize the symmetry of the diffraction 
rings in both the vertical and horizontal directions. Af- 
ter alignment, the membrane was exactly perpendicular 
to the incident beam and the long axes of the pores were 
parallel to the incident beam. This defined 9s = 90° be- 
tween the incident beam and the short axis of the pores. 
The membrane was then rotated to 0= = 10 ± 1° from 




Pores 



CCD detector 



FIG. 3: Geometry for SAXS measurements: X-rays were in- 
cident at 6a — 10° from normal to the nanopore long axis. 
Scattered intensity was measured with a fixed position CCD 
detector down-stream. The m and q describe the coordinate 
space and wavevector transfer of the nanopore, respectively. 
Note that the geometry shown here has been rotated about 
the incident beam by 90° for clarity. 



the incident beam, i.e. the pore long axis was then at 
80° to the incident x-rays (see Fig.|31). This geometry al- 
lowed a bulk measurement of the nanoparticle/nanopore 
system while maximizing the wavevector transfer along 
the nanopore long axis, Qz, that was probed. 

It is important to note that the geometry shown in 
Fig. O is rotated by 90° about the incident beam for ease 
of viewing. In the figure, the long axis of the pores is 
depicted as approximately perpendicular to the incident 
beam and vertical with rotations 9^ in the vertical plane. 
In the actual experiment, the only difference was that the 
long axis of the pore was horizontal (though still perpen- 
dicular to the incident beam) and rotations 9s were done 
in the horizontal plane. 



IV. VOLUME OF ADSORBED SOLVENT 

A measure of absorbed volume of solvent in the pores 
was calculated from the small angle powder diffraction 
peaks mentioned above. As liquid was condensed into the 
pores, the electron density contrast in the pores decreased 
relative to that of the dry pores, thus reducing the scat- 
tering intensity of the diffraction peaks, Ipeak{AT), rela- 
tive to the dry peaks, Ipeak{dry). Neglecting absorption 
corrections, the added solvent volume, Vuq, was then 
approximately determined from the lowest order peak 
(small angle approximation) viaSi: 

Vii,(AT) cx ^J Ipeak (dry) - ^ Ipeak (AT) (1) 

Plots of the added volume derived from the (10) diffrac- 
tion peak are shown in Fig. ^ The data represent three 
different thermal cycles: two cooling cycles (large AT 
to small AT) and one heating cycle. The data have 
been normalized to the saturated (pores completely filled 
with liquid) volume for each curve. There are two main 
features in this figure: the first is the presence of a 
sharp transition (from low relative volume to saturation 
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FIG. 4: (Main) Volume adsorption(A, □) / desorption(O) 
curves as a function of AT. Dashed lines are added as a 
guide for the eye. Note the strong hysteresis upon cycling. 
Vertical dotted lines indicate location of capillary transition 
for empty (no nanoparticles) pores predicted by the Kelvin 
equation, A/i ~ y/pRc- Data is normalized to the liquid sat- 
uration volume of the pores. (Inset) Plot showing the powder 
diffraction peaks from the 2D nanopore packing. Volume in- 
formation was extracted from the first peak (the (10)). 



and reverse) for both adsorption and desorption curves 
and second, there is a marked hysteresis which is repro- 
ducible (at least for the cooling cycle). For both sets 
of curves, the nanopore/nanoparticle system is saturated 
with toluene for AT < 4 K and has negligible amounts of 
toluene for AT > 30 K. These two regions will be referred 
to throughout the rest of the paper as the "saturated" 
and "dry" regions, respectively. 

The simplest theory to describe the observed capillary 
transition is the Kelvin equation: 



adsorption 
orption 
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where Rc is the cylindrical nanopore radius, 7 — 
28.4 mN/m is the surface tension of bulk toluene, and 
ni — p/ra^t = 9 x 10^^ moles/cm'^ (mass density di- 
vided by molecular weight) is the molar density of the 
toluene. The chemical potential at the capillary transi- 
tion for the adsorption and desorption differ by a factor 
of two due to different physical mechanisms for adsorp- 
tion and desorption. Pores are expected to fill by coax- 
ial film growth (curvature of l/Rc) and empty from the 
ends of the pores via a spherical meniscus (curvature of 
2/Rc) that travels the length of the poresS? This argu- 
ment provides a motivation for the presence of the hys- 
teresis, even though the system studied is clearly more 
complicated due to the presence of the nanoparticles. Ac- 
cording to the Kelvin equation, the observed capillary 
filling transition in Fig. 01 would indicate an effective ra- 
dius for the pores of about 4 nm, almost a factor of four 
smaller than the actual pore radius of 15 nm. For simi- 



lar nanoporous alumina samples with no nanoparticles, it 
was found that the Kelvin equation provides a good pre- 
diction of the capillary transitions^. Significantly higher 
values of AT for the capillary transition indicated that 
there was a real reduction in the effective radius of the 
pores due to the presence of the nanoparticles. Assum- 
ing a monolayer thickness due to the nanoparticles of 
2.4 nm for the Au core plus 1.2 nmx2 for the organic 
OT shell22«:! 5 nm, this would bring the effective radius 
down to about 10 nm. Additionally, roughness from the 
nanoparticle monolayer surface contributes to the tran- 
sition shift by increasing the amount of adsorbed liquid 
at large AT^. 



V. NANOPARTICLE-NANOPARTICLE 
SCATTERING 

A. Elliptical Transforms 

In addition to the low angle powder diffraction from 
the 2D hexagonal pore packing, a ring-like structure was 
observed at larger angles [q ~ 0.18 , see Fig. EJ, 
which corresponds to nanoparticle interference scattering 
and contains information about the local packing struc- 
ture. This scattering ring underwent dramatic changes 
with the gradual addition, and subsequent removal, of 
liquid solvent in the pores (see Fig.EJ. No sharp diffrac- 
tion spots were observed upon rotation of the sample 
through 90° about the short axis of the pore, indicating 
that the nanoparticle packing must be powder-like with 
only short-range order. 

To interpret these particle-particle scattering results, 
it was necessary to transform the scattering intensity 
from lab (CCD detector) coordinates into the coordi- 
nates relative to the nanopore axis as shown in Fig. 
The transform is given by the constraints of the scat- 
tering geometry (see Fig. The wavevector transfer, 

Q = [Qx,Qy,Qz] is: 



Q — k'inc ^scatt 

k 



x,L — \/ x-^ -\- z-^ + L^, z 



-[a;,0,z]; T>r 



(4) 
(5) 

(6) 



Transforming to coordinates of the nanopore (here we 
use (f, lowercase for clarity), for 6s ^ 7r/2: 



q ■ rriz kz cos 9s 



qy 



\m 

q- my 



2 L 
kz sin 8s 



y\ 



q ■ rux 
k\/ 



kx 

T 



z^sin^^ 



(7) 
(8) 
(9) 

-; T > r; z » Acot6's (10) 
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FIG. 5: (Color) (Left) Schematic of elliptical transforms of 
scattering intensity data. Lines of constant q,. are ellipses in 
CCD coordinates where the scattering (CCD) plane intersects 
the cylindrical surfaces of constant qr- (Right Top) CCD im- 
age of scattering intensity at AT = 30 K. The direct beam is 
in the lower left corner. The triangular region here is due to 
an attenuator used to block the intense pore-pore scattering 
at low q. (Right Bottom) The same data set transformed to 
the q coordinates relative to the nanopore axis. 



Here, x — [x,y,z\ {\\x\\ — r) and m = [mx.,my,mz\ are 
the cartesian coordinates of the CCD detector (lab) and 
nanopore, respectively, and A — \/ x'^ + + L'^. The 
wavevector transfer in the pore coordinates is denoted q, 
with magnitude k = 2t:/X. In the above approximation, 
the intersection of the cylindrical surfaces of constant qr, 
and the CCD plane, are ellipses. Fig. \E\ shows typical 
scattering intensity data (upper right) from AT = 30 K 
(dry pores). The scattering intensity is transformed into 
the qz,qr coordinates of the nanopore (lower right). The 
transform excludes data up to 0s = 10° from the qz axis 
(aligned along the nanopore long axis) as a result of the 
incident angle of the x-rays to the pore axis. Images 
were taken with the detector off-center to maximize the 
recorded g-range. The triangular region at low q is a 
result of an attenuator necessary to reduce the intense 
scattering of the powder diffraction peaks associated with 
the hexagonal pore-pore packing. The intense scattering 
concentrated along the qr axis was mainly a result of scat- 
tering described by the pore form factor which is narrow 
in qz due to the length of the pores. 

The intensity distribution along the nanoparticle scat- 
tering ring changed significantly with AT. For the dry 
pores, AT = 30 K (Fig. El upper left) a sharp, strongly 
asymmetric ring was present, brighter near the qz axis 
than the qr axis. Here the data have been tiled to all four 
quadrants to simulate the full ring for ease of viewing. 
The strong feature along the qr axis was due to the sum 
of the scattering of the individual pores. The ring asym- 
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FIG. 6: (Color) Scattering intensity distributions at four rep- 
resentative AT. Data (intensity log scale) were background 
subtracted, then tiled into all quadrants to simulate a full ring. 
(Top Left) Data from AT = 30 K (dry): the scattering in- 
tensity is an asymmetric ring structure from particle-particle 
scattering. (Top Right) Data from AT = 8 K (adsorption 
curve): note the ring is broader with a smaller radius. (Bot- 
tom Right) Data from AT = 6 K (saturated with liquid): 
here the ring is even broader with an even smaller radius. 
(Bottom Left) Data from AT = 12.4 K (desorption curve): 
note the presence of an isotropic scattering ring just inside 
the asymmetric ring. 



metry indicates a structure that is preferentially aligned 
along the nanopore axis. A model to explain this scatter- 
ing will be presented in the next section. As liquid was 
added to the pores (adsorption curve) there was a gen- 
eral trend that the ring became broader and its radius 
decreased (see Fig. El AT = 30 K ^ 8 K ^ 6 K). This 
was an indication that the nanoparticle nearest neigh- 
bor spacing had increased and the ordering was reduced 
in comparison with that of the dry case. Upon grad- 
ual removal of the liquid (desorption curve) there was a 
qualitatively different feature that appeared (see Fig. El 
lower left, AT 12.4 K). For this value of AT, the in- 
ner portion of the asymmetric ring was nearly isotropic. 
This isotropic scattering was an indication of an addi- 
tional structure of nanoparticle aggregates or clusters. 
With removal of nearly all of the liquid achieved at the 
highest AT, the scattering again showed a sharp, well 
defined, asymmetric ring (almost identical to Fig. El up- 
per left, AT = 30 K), indicating that the self-assembly 
was reversible. The formation of the clusters indicated 
that there was a form of hysteresis in the self assembly 
process as well as in liquid adsorption and desorption. 
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B. Nanoparticle Tiling Model 

We propose in this section a simple model to describe 
the scattering from monolayers of Au nanoparticles along 
the walls of the cylindrical nanopores. Motivation for 
a self-assembled cylindrical monolayer comes from anal- 
ogy with self-assembled nanoparticle monolayers on flat 
substratesi and the attractive van der Waals (vdW) 
forces between the particles and between the particles 
and the walls of the nanopore. In the Derjaguin ap- 
proximation, where the separation, Z?, is much less than 
the nanoparticle Au-core radius, Rg, the vdW forces be- 
tween neighboring nanoparticles, Fnp-np, and between 
a nanoparticle and the alumina pore wall, Fnp-waii, 
are»ia*S. 



F, 



np—wall 



121)2 

6D^ 



(11) 

(12) 



Here, A„„ and Anw sue the Hamaker constants for the 
nanoparticle-nanoparticle interaction and nanoparticle- 
nanopore wall respectively. The Hamaker constants are 
reduced in the presence of a mediating solvent, but the 
forces remain attractive. Steric repulsions, due to the 
OT coating, prevent irreversible nanoparticle aggrega- 
tion. The attractive vdW forces should partially confine 
the particles to remain in a monolayer near the walls of 
the pores. Additionally, inter-particle ordering should be 
reduced by pore wall roughness and nanoparticle polydis- 
persity.^^a The lack of strong diffraction peaks confirmed 
that the nanoparticle order was short-range only. There 
were two distinct length scales involved, the nanoparticle 
diameter and the nanopore diameter, that differ by an or- 
der of magnitude (~ 3 nm to ~ 30 nm). Thus, there were 
two independent scattering regimes: one for low angles, 

qRc ~ 1, and another for higher angles, qRg ~ 1, where 
Rs is the average nanoparticle core radius and Rc is the 
average nanopore radius. For the low angle scattering, 
which has already been discussed above in the context of 
added liquid volume, the scattering due to interference 
from neighboring nanoparticles was ignored for this rea- 
son. For the scattering at higher angles, the structure 
factor associated with the nanopore packing should have 
decayed to approximately unity (since any pore-pore or- 
dering must be short range in 2D). This high q regime 
may then be written as the sum of scattering from the 
pores (as in the low q region, but where now the struc- 
ture factor is unity) plus a term that describes the local 
particle-particle scattering. The scattered intensity in 
the two regimes is given by: 



I{q) « \Sp{q)iF,iq) + Fp{q))\-, qRc 
liq) « \Fpiq) + Snpiq)Fnp{q)\^; qRs 



1 



(13) 
(14) 
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FIG. 7: (Left) Schematic of the tiling of the cylindrical 
nanoparticle monolayer. The monolayer is broken up into 
small regions, or tiles, which are approximately flat. The 
scattering is then powder averaged over all tiles at a given 
orientation (j) from all pores. The total scattering is a sum 
over all such orientations, (j) £ [0, Stt]. The azimuthal angle is 
^ (not to be confused with the sample rotation angle 9s). 



respectively. Fp, Fs, and Fnp are the form factors associ- 
ated with the scattering from the nanopores, an average 
nanoparticle monolayer, and individual nanoparticles re- 
spectively. Since the nanopores are about four orders 
of magnitude longer than their diameter, the scattering 
associated with them will occur at small q^. With the 
pores almost perpendicular to the direct beam, the in- 
tersection of this scattering with the image plane will 
be such that scattering associated with the pores will be 
only along the qr axis. The scattering associated with 
the form factor of the pores is complicated by polydis- 
persity, roughness and the presence of liquid, thus it is 
treated only as background scattering for this analysis. 
Strong scattering away from the qr axis must then be 
attributed to nanoparticle-nanoparticle scattering. Thus 
away from the g,. axis, we may treat the 2nd term in 
equation (|14|l as dominant and the first term (and cross 
term) as a background. 

The nanoparticle form factor is described by: 



np I 



G(g,i?„Afl). 



(15) 



where G{q, Rs, An) is given by the Shulz distribution^ 
for polydisperse spheres with radius Rs and Ar defined 
as the root mean square deviation. For qRg ^ 1, the 
nanoparticle form factor is given by the Guinier Law2^. 



Giq) « exp(-g2^, V5) 



(16) 



Here, Sp and Snp are the structure factors associated with 
the local packing of the nanopores and the nanoparticles 



A "tiling" model for calculating Snp (q) , valid for qRg ~ 
1 is explored below. To model the nanoparticle structure 
factor for high q (where nanoparticle-nanoparticle inter- 
ference scattering dominates) , the monolayer of particles 
was broken up into small regions, or tiles, that were ap- 
proximately flat, as shown in Fig. [T] The approximation 
was made that the scattering between individual tiles was 
independent and that the radius of curvature of the pore 
is large compared with the inter-particle distance. The 
first approximation is valid when the dimensions of the 



tiles are on the order of the average correlation length 
(order parameter), ^, within the tile or smaller. 

The scattering was then powder-averaged over all tiles 
of a given orientation 0. Thus the scattering intensity 
from a particular set of tiles at a given was approxi- 
mately the same as the powder-averaged scattering from 
a flat monolayer of nanoparticles with an order parame- 
ter ^, and nanoparticle nearest neighbor separation (i„„. 
A (mathematically) simple Lorentzian model was used to 
approximate the structure factor, S**^, for each of these 
orientations: 

|5t((?,0,V')l'« (17) 

M 

1 + {q\/ (cos('i/') - cos(0) sin(V')) - 2n/ d„„^ 

Here Iq is an adjustable scale parameter with Ig/dnn be- 
ing proportional to the number of scatters and ip is the 
azimuthal angle relative to the long axis of the nanopore 
(not to be confused with the sample rotation, 9s). Note 
that if ^ > dnn then ^ may be thought of as the correla- 
tion length within the tiles, along the walls of the pore. 
For ^ < dnn then the particles are better described as a 
2D dilute, gas-like phase. This analysis only considered 
the first order peak due to the q-range of this experiment, 
but higher orders would be treated similarly. The total 
structure factor, Snp, can be calculated by integrating 
over all (f>. 

\Snp{q,0)\^ = r \St^{q,cl,,0)fd^ (18) 
Jo 

The above structure factor, along with the nanoparti- 
cle form factor given by the Shulz distribution, was used 
to fit the high q data of the Au-core nanoparticles in 
the nanoporous alumina. The CCD data were first flat- 
field corrected, then background (taken with the sample 
rotated out of the beam) was subtracted to account for 
scattering from our environmental chamber windows plus 
scattering from the toluene vapor. The data were then 
transformed into the qr,qz coordinates of the nanopore 
as described above. All of the data sets (at different 
AT) fit well with the above described tile model with 
a monotonic decaying background (which included scat- 
tering from both the porous membrane and the liquid) 
with the exception of the data at AT = 12.4 K (desorp- 
tion curve) where an additional isotropic scattering term 
was necessary for good agreement. This cross-section of 
the isotropic scattering ring was fit using a Lorentzian 
with similar parameters as for the monolayer {^,Io, and 
dnn)- For this model, the number of scatters should be 
proportional to Iq / d^n ■ 

For each data set of a given AT, a series of represen- 
tative radial slices was chosen at ■(/'„ = 11 x n° from the 
qz axis (n an integer from one to seven), and intensity 
as a function of q was graphed for each slice (see Fig.lSJ. 
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FIG. 8; (Color) (Right Top) Schematic of representative ra- 
dial slices (dashed lines) taken for fitting with the tile model 
with ?/>„ = 11 X n°, n an integer from 1 to 7. (Left) Data 
(□) at AT = 30 K from each slice (color coded with dashed 
lines) along with fits (solid lines) from the tile model plus 
independent monotonically decaying, positive backgrounds. 
Error bars are smaller than data symbols. (Bottom Right) 
Data (all at = 11°) from four different AT: □ = 30 K, 
A = 8 K (absorption), = 4.4 K(saturated), and Q — 
12.4 K(desorption); solid lines are fits. AT — 12.4 K is qual- 
itatively different and requires an isotropic component in ad- 
dition to the tile model to fit properly. 



Each radial slice was fitted independently with the tiling 
model plus a monotonically decaying background. As 
stated above, the main contribution to the background 
was from the nanopore membrane which increased as ipn 
approached 90°. This was seen in the fits and thus a slice 
at 88° was not used due to the fact that the background 
was greater than the particle-particle scattering. After 
fitting all of the slices independently, average parameters 
and uncertainties for each AT were calculated from the 
mean and standard deviation of each of the three fit pa- 
rameters: /q, and dnn- Since the uncertainties were 
small in comparison to the average value of the parame- 
ter, the average parameters were treated as "best fit" or 
global parameters for each AT. 

The three "best fit" parameters from the tiling model 
as a function of AT are plotted with uncertainties in 
Fig. El Also plotted in each of the three plots are the pa- 
rameters from the isotropic component at AT = 12.4 K 
(■). The AT axis is read from left to right; to the left of 
zero are adsorption (cooling) data, to the right of zero is 
desorption (heating) data. 

The data indicated that the number of nanoparti- 
cles in the monolayer, probed by the scattering vol- 
ume, remained constant as liquid was added and sub- 
sequently removed (see Fig.|31 top) with the exception of 
AT — 12.4 K where an additional isotropic component 
was necessary. For this data point, the number of par- 
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FIG. 9: Plots of the "best fit" tile model (Q) and isotropic 
(■) parameters and uncertainties as a function of AT (data 
to the left of zero are from adsorption curves, to the right are 
from desorption curves). The physical parameters are (top 
to bottom): number of scatters oc Iq /<!„„, nearest neighbor 



spacmg 



and order parameter 5. Lines are added as a 



clusters would dissolve and reassemble into a monolayer 
structure on the nanopore wall. This is indicated in the 
data by the disappearance of the isotropic scattering as 
AT is increased. 

The plots in Fig. [51 middle and bottom, demonstrate 
that with the addition of liquid there was a shift to larger 
nearest neighbor separation distance with a maximum of 
4.5 nm (compared to 3.6 nm for dry), accompanied by 
a decrease in the order (in the plane of the cylindrical 
monolayer). Upon removal of the liquid solvent, the par- 
ticle separation decreased again to the dry value, and the 
ordering increased back to almost the initial dry value. 
For relative added liquid solvent volume amounts up to 
about 0.5 (normalized to saturation; see both Fig. 2|and 
Fig. O the order parameter was about twice that of the 
nearest neighbor spacing and may be interpreted as an 
in-plane correlation length. For higher relative amounts 
of solvent, Vuq >0.5, the order parameter was less than 
the nearest neighbor spacing, indicating a 2D dilute gas- 
like phase. 

The value of 3.6 nm for the nearest neighbor spacing 
for the dry system (large AT = 30 K) was less than 
2{Rs + tor) ~ 4.8 nm, where tor ~ 1-2 nm is the thick- 
ness of the OT shell, indicating interdigitation of the shell 
ligands of neighboring particles. As liquid was absorbed 
into the organic shells (with decreasing AT) the vdW 



attraction between nanoparticles would decreas( 



.12.13 



•el- 



guide to the eye. 



ative to the dry system, while the osmotic pressure^ 
between the ligand chains increased, driving the parti- 
cles apart. This separation increase may have been fur- 
ther facilitated by gaps or voids in monolayer coverage. 
As liquid filled these gaps, the nanoparticles (undergoing 
thermal motion) could move into this new volume. Re- 
moval of the liquid by increasing AT reduced the repul- 
sive osmotic pressure between the particles and increased 
the attractive vdW forces, reducing the particle-particle 
separation. 



VI. CONCLUSIONS 



tides in the monolayer decreased. It is likely that some 
of the missing particles went into the formation of the 
isotropic clusters. The fact that this data point is higher 
than the average monolayer number might be indicative 
that the clusters fill the volume more effectively than the 
particles in the monolayer. 

The formation of these clusters is a surprising phe- 
nomenon. Additional isotropic scattering was seen only 
during desorption and then disappeared upon removal 
of additional liquid. The formation of clusters, or ag- 
gregates may have been linked to the progression of a 
spherical meniscus through the pores, starting from the 
endsi^ As the meniscus travelled through the pore, it 
would tend to drag particles with it. Another possibil- 
ity would be the formation of bubbles within the liquid. 
Upon complete removal of liquid from the system these 



Low angle measurements established the relative 
amount of liquid in the nanopores as a function of AT. 
A capillary filling transition occurred between 8 K and 
12 K, which was about four times less saturation than 
what is expected via the Kelvin equation for toluene ab- 
sorption in pores with no nanoparticles. In addition to 
the shift, marked hysteresis was also observed. 

Fits to the data indicated that the nanoparticles as- 
sembled in a cylindrical monolayer aligned along the 
pore. From the geometry of this structure we conclude 
that the particles were near the walls, which was physi- 
cally sensible due to the vdW attraction of the nanopar- 
ticles to the nanopore wall. As liquid was added by re- 
ducing the relative chemical potential A/i, most of the 
particles remained in the monolayer structure. Also, the 
nearest neighbor separation distance increased and the 
correlation length within the monolayer decreased. The 
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process was reversible: upon removal of the liquid, the 
nanoparticle nearest neighbor distance decreased to the 
initial dry value and the ordering increased to almost the 
dry value. 

In addition to the cylindrical monolayer, there was 
evidence of the formation of isotropic clusters during 
desorption of the liquid solvent. This phenomenon could 
have been related to the process by which the liquid 
emptied from the pore, namely that it was likely to occur 
from the ends of the pores. This might also explain why 
this structure was not seen during the adsorption cycle. 
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